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ABSTRACT 



Chandra has recently observed IE 0657-56, a hot merg- 
ing system at z = 0.3 (the "bullet" cluster), for 500 ks. I 
present some of the findings from this dataset. The clus- 
ter exhibits a prominent bow shock with M = 3.0 + 0.4 
(one of only two known M » 1 shock fronts), which we 
use for a first test of the electron-ion equilibrium in an 
intergalactic plasma. The temperatures across the shock 
are consistent with instant shock-heating of the electrons; 
at 95% confidence, the equilibration timescale is much 
shorter than the collisional Spitzer value. Global prop- 
erties of IE 0657-56 are also remarkable. Despite being 
extremely unrelaxed, the cluster fits well on the Lx~T re- 
lation, yet its total mass estimated from the M-T relation 
is more than twice the value measured from lensing. This 
is consistent with simulations predicting that in the mid- 
dle of a merger, global temperature and X-ray luminosity 
may be temporarily boosted by a large factor. 

Key words: galaxies: clusters: individual (1E0657-56) 
— plasmas — X-rays: galaxies: clusters. 



1. INTRODUCTION 



IE 0657-56 is the most interesting cluster in the Uni- 
verse, as officially confirmed by Chandra Peer Review 
(Anonymous 2003). This system, located at z = 0.3, 
has the highest X-ray luminosity and temperature and 
the most luminous radio halo of all known clusters. It 
is also a spectacular merger occurring almost exactly in 
the plane of the sky (Markevitch et al. 2002), and con- 
tains one of only two known cluster shock fronts with a 
Mach number substantially greater than 1 (the other one 
is A520 with M = 2). Chandra has recently observed 
it with a 500 ks total exposure. The image from that 
dataset is shown in Fig. 1. It shows a prominent bow 
shock preceding a small, cool "bullet" subcluster flying 
west after passing through a core of a bigger cluster and 
disrupting it. In this paper, two interesting (preliminary) 
findings from this new dataset are presented, one regard- 
ing the global properties of the cluster, and another based 
on the high-resolution electron temperature profile across 



the shock front. All errors are 68%; the assumed cosmol- 
ogy is h = 0.7, Qo = 0.3 and Q A = 0.7. 



2. AN OVERHEATED CLUSTER 



A single-temperature fit to the overall ACIS spectrum for 
the cluster (excluding a small region around the bullet 
for consistency with the "cooling flow corrected" tem- 
peratures for nearby clusters) is T = 14.1 ± 0.2 keV. As 
shown in Fig. 2, it fits perfectly on the Lx - T relation 
for local clusters (Markevitch 1998) after a correction for 
its redshift evolution (Vikhlinin et al. 2002). However, 
its total mass estimated from an X-ray M500 - T rela- 
tion (Vikhlinin et al. 2005b; Kotov & Vikhlinin 2005), 
M500 = 1.9 x 10 I5 M o , is a factor of 2.4 higher than the 
value within the same radius estimated from weak lensing 
(Clowe et al. 2004). Given the ongoing violent merger, 
this is not unexpected — simulations have predicted that 
during a major merger, the cluster may experience a tran- 
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Figure 2. IE 0657-56, shown as open square, overlaid 
on the local Lx — T relation ( with its best-fit power law 
shown by dashed line), after a correction of Lx of this 
cluster for the redshift evolution of the relation. 



Figure 1. 500 ks Chandra ACIS-I image of IE 0657-56 in the 0.8-4 keV band. 



sient boost of temperature by a factor of several, lasting of 
order 0.1 Gyr around the moment of the subcluster core 
passage (Randall et al. 2002; Rowley et al. 2004). In the 
course of this rapid change, the cluster moves approxi- 
mately along the L-T relation. We know from the shock 
velocity (Markevitch et al. 2002) that the core passage in 
IE 0657-56 has indeed occurred about 0.15 Gyr ago, so 
this cluster appears to illustrate precisely this short-lived 
phenomenon. 

It is interesting to try and recover a "pre-merger" tem- 
perature of the main subcluster. The gas before (west of) 
the shock front should not yet know anything about the 
merger, and stay undisturbed in the gravitational potential 



of the main cluster. We have therefore divided the cluster 
into the post-shock and pre-shock regions (approximately 
along the bow shock) and extracted radial temperature 
profiles in these regions, both centered at the mass cen- 
troid of the main subcluster (Clowe et al. 2004). The re- 
sulting profiles are shown in Fig. 3. If before the merger, 
the main cluster had a declining radial profile similar to 
that in most clusters (Vikhlinin et al. 2005a), extrapolat- 
ing the "pre-shock" profile inwards would give an aver- 
age temperature around 10 keV or less, compared to the 
present 14 keV. For such a temperature, M500 from the 
M - T relation is within a factor of 1.5 of the lensing 
mass, which is within the uncertainty of the Clowe et al. 
(2004) lensing measurement. (The lensing mass would 
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Figure 3. Radial temperature profiles for the disturbed 
region behind the shock ( crosses without bars) and for the 
undisturbed region in front of the shock (barred crosses). 
Both profiles are centered on the main subcluster. The 
vertical dashed line approximately shows the shock front; 
the horizontal line is the average temperature. 

imply T - 8 keV.) 

The "overheated" cluster IE 0657-56 may thus serve as 
a cautionary example for projects involving X-ray sur- 
veys of distant clusters, where total masses are proposed 
to be derived from low-statistics X-ray data using tem- 
perature profiles, average temperatures or even fluxes. If 
IE 0657-56 was placed at z — 1 and observed with Chan- 
dra with a relatively long 100 ks exposure, its extremely 
unrelaxed state would be difficult to detect from the X- 
ray data (image, Lx - T relation, etc.), especially if the 
merger was not oriented so fortunately in the plane of the 
sky. Its mass estimate would of course be significantly 
wrong. 



3. ELECTRON-ION EQUILIBRIUM 

The bow shock in IE 0657-56 offers a unique experimen- 
tal setup to determine whether electrons in the intracluster 
plasma are directly heated by shocks, or compressed adi- 
abatically and then heated to an equilibrium temperature 
via collisions with protons (that are heated dissipatively 
by the shock). The collisional equilibration occurs on a 
Spitzer timescale (e.g., Zeldovich & Raizer 1966) 

r ep =2xl0 8 yr( ^—Y^Jf-f 2 (1) 

p 3 \10- 3 cm- 3 / \10 8 K/ 

(note that, conservatively for our measurement, this is 3 
times shorter than the formula given in Sarazin 1988). 
We cannot measure T-, in X-rays, only T e . However, be- 
cause the shock in IE 0657-56 propagates in the plane 



Figure 4. X-ray brightness profile across the shock front. 
The line shows the best-fit model (a projected sharp 
spherical density discontinuity at the shock). 

of the sky, we can accurately measure the gas density 
jump across the front and use it to predict the post- 
shock adiabatic and shock-heated (dissipative) electron 
temperatures from the pre-shock temperature (using the 
adiabat and the Rankine-Hugoniot jump conditions, re- 
spectively), and compare it with the observation. Fur- 
thermore, we also know the downstream velocity of the 
shocked gas flowing away from the shock. This flow ef- 
fectively unrolls the time dependence of the electron tem- 
perature along the spatial coordinate for us. The Mach 
number of the shock is conveniently high, so that the 
adiabatic and dissipative electron temperatures are suf- 
ficiently different for us to distinguish (e.g., for M < 2, 
they would be almost the same). It is also not a strong 
shock, for which the density jump would just be a fac- 
tor of 4 (for ideal monoatomic gas) and would not let 
us directly determine M. Furthermore, the distance trav- 
eled by the post-shock gas during the time given by eq. 1, 
Ax 230 kpc = 50", is well-resolved by Chandra. The 
statistical quality of the 500 ks dataset is just sufficient — 
in fact, this test was the main science driver for this long 
observation. 

Fig. 4 shows a 0.8-4 keV surface brightness profile in 
a narrow sector across the shock, centered on the center 
of curvature of the front. The inner bump is the bullet 
(its boundary is a "cold front"). The edge at 90" is the 
shock front. There is also a subtle secondary edge be- 
tween these main features, which may also be seen in the 
deep image. It is unrelated to the shock, so we should 
take care to exclude it from any fits aimed at deprojecting 
the shock temperatures and densities. The line in Fig. 4 
shows a best-fit model consisting of the projected abrupt 
spherical density jump at the shock (by a factor of 3.0), a 
power-law profile inside the shock and a beta-model out- 
side. The fit is perfect (the inclusion of the secondary 
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Figure 5. Left: projected temperature profile in a narrow sector across the shock. Two crosses in the shock region with 
lower temperatures correspond to an apparent additional edge-like structure and are not used. Vertical lines show the 
boundaries of the cool bullet and the shock; dashed line shows the average pre-shock temperature. Right: deprojected 
temperatures for the two post-shock bins overlaid on the model predictions (with error bands) for shock-heated and 
adiabatic electron temperatures. The velocity shown is for the post-shock gas relative to the shock. 



edge does not affect it). From this density jump, we de- 
termine M = 3.0+0.4, which corresponds to a shock (and 
bullet) velocity of 4700 km s" 1 . 

In Fig. 5 (left), we show a projected temperature profile 
in the same sector across the shock. There is a clear jump 
of the electron temperature. At the secondary edge men- 
tioned above, the temperature goes down, which probably 
indicates residual cooler gas from the subcluster located 
ahead of the bullet. Therefore, we can only use the two 
bins closest to the front. Because we know the gas den- 
sity profile across the shock, we can accurately subtract 
the contributions of the cooler pre-shock gas projected 
into the post-shock bins, assuming spherical symmetry. 
The large brightness contrast at the shock helps make this 
subtraction robust. 

The deprojected values are shown in Fig. 5 (right). They 
are overlaid on the two models (gray bands), one as- 
suming instant equilibration (i.e., electrons are heated at 
the shock), and another assuming adiabatic compression 
and subsequent equilibration with protons on a timescale 
given by eq. 1 . (The plot assumes a constant post-shock 
gas velocity, which of course is not correct, but we are 
interested only in the immediate shock vicinity.) The de- 
projected gas temperatures are so high for Chandra that 
only their lower limits are meaningful. The temperatures 
are consistent with instant heating; the "adiabatic" model 
with the Spitzer timescale is excluded at a 95% confi- 
dence. 

A few sanity checks have been performed. The high tem- 
peratures are not an artifact of the deprojection, because 
the projected temperatures in those two bins are already 



higher than the models. We also considered the possibil- 
ity of a non-thermal contamination. The cluster is known 
to possess a radio halo (Liang et al. 2000) which has an 
edge right at the shock front (Markevitch et al. 2002). 
Therefore, there may be an inverse Compton contribu- 
tion from relativistic electrons accelerated at the shock. 
However, the power-law spectrum of such emission for 
any M would be softer than thermal, so should not bias 
our measurements high. We also extracted a tempera- 
ture profile in another sector of the shock (away from the 
nose), where M 2, and made a similar comparison. The 
recovered post-shock temperatures are lower than those 
in Fig. 5 and again in agreement with the model predic- 
tions (either of the two models; for such M, the differ- 
ence between them is small). Thus, albeit at a relatively 
low significance (95%), we conclude that the electron-ion 
equilibration should be much faster than collisional. 

It is of course unfortunate that a cluster with such a per- 
fect geometric setup and Mach number is so hot that 
Chandra can barely measure the post-shock tempera- 
tures. However, there is no choice of shock fronts and 
improvements upon the above measurement are unlikely 
in the near future. 

There hasn't been a measurement of the electron-ion 
equilibration timescale in the intergalactic medium be- 
fore. In the solar wind plasma, the equilibration is be- 
lieved to be fast compared to the collisional timescale. 
For supernova remnants, which have strong shocks, con- 
clusions vary between different objects (e.g., Rakowski 
2005). Plasma interactions have been suggested as a fast 
equilibration mechanism for the solar wind, and it proba- 
bly applies for the intracluster plasma as well. 



4. SUMMARY 



From the extra-long 500 ks Chandra observation of 
IE 0657-56, we found that this cluster is observed at a 
very special, short-lived stage, when its temperature and 
luminosity are temporarily boosted by the merger by sig- 
nificant factors. The total cluster mass estimate from the 
X-ray M -T relation turns out to be more than two times 
higher than the (presumably) true mass determined by 
lensing, indicating a very strong deviation from hydro- 
static equilibrium. If this cluster were at high z and not so 
well exposed, it would be difficult to detect its disturbed 
state, thus it is a cautionary example for future high-z sur- 
veys. 

The temperature profile across the shock offers a first 
test of the electron-ion equilibrium in an intracluster 
plasma. The temperatures indicate that electrons are in- 
deed quickly heated at the shock. The slow collisional 
(Spitzer) equilibration rate is excluded at the 95% confi- 
dence. 

This unique cluster, well-exposed by Chandra, is the sub- 
ject of several other ongoing studies, the results of which 
will be reported soon. 



REFERENCES 

Clowe, D., Gonzalez, A., & Markevitch, M. 2004, ApJ, 
604, 596 

Kotov, O., & Vikhlinin, A. 2005, astro-ph/051 1044 

Liang, H., Hunstead, R. W., Birkinshaw, M., & An- 
dreani, P. 2000, ApJ, 544, 686 

Markevitch, M. 1998, ApJ, 504, 27 

Markevitch, M., Gonzalez, A. H., David, L., Vikhlinin, 
A., Murray, S., Forman, W., Jones, C, & Tucker, W. 
2002, ApJ, 567, L27 

Rakowski, C. E. 2005, Adv. Space Research, 35, 1017 

Randall, S. W., Sarazin, C. L., & Ricker, P. M. 2002, 
ApJ, 577, 579 

Rowley, D. R., Thomas, P. A., & Kay, S. T. 2004, MN- 
RAS, 352, 508 

Sarazin, C. L. 1988, X-ray emission from clusters of 
galaxies (Cambridge University Press) 

Vikhlinin, A., VanSpeybroeck, L., Markevitch, M., For- 
man, W. R., & Grego, L. 2002, ApJ, 578, LI 07 

Vikhlinin, A., Markevitch, M., Murray, S. S., Jones, C, 
Forman, W., & VanSpeybroeck, L. 2005a, ApJ, 628, 655 

Vikhlinin, A., Kravtsov, A., Forman, W., Jones, C, 
Markevitch, M., Murray, S. S., & Van Speybroeck, L. 
2005b, astro-ph/0507092 

Zeldovich, Y. B., & Raizer, Y. P. 1967, Physics of Shock 
Waves (New York: Academic Press) 



